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EXTENDING  THE  CAPABILITIES  OF  THE  NATICK  TERRAIN 
ANALYSIS  SYSTEM 


1.  Introduction. 

The  Natick  Terrain  Analysis  System1"3  was  designed  to  satisfy  the  need 
for  a  more  scientific  method  of  designing  camouflage  patterns  and  coloration 
based  on  actual  terrain  reflectance  data.  This  report  describes  the  three- 
phase  expansion  of  the  Terrain  Analysis  System  to  include  reflectance  data  in 
the  near-infrared,  to  allow  for  the  analysis  of  video  data  obtained  through  an 
image  intensifier  or  thermal  imager,  and  for  the  interactive  creation  and 
evaluation  of  camouflage  patterns. 

The  Phase  I  effort  can  be  subdivided  into  six  tasks: 

a.  Modify  the  Video  Data  Processing  software  to  allow  for  wavelengths  in 
the  range  of  720  nanometers  (nm)  to  900  nm. 

b.  Convert  the  400-900  nm  data  to  L*  values  using  published  photocathode 
sensitivities  and  nighttime  illumination  data. 

c.  Plot  the  reflectance  curve  of  a  single  pixel  representative  of  each 
clustered  domain. 

d.  Produce  varying  density  symbol  plots  of  the  clustered  domains. 

e.  Model  three  dye,  false  color,  infrared  film. 

f.  Plot  the  clustered  domains  on  the  LA-100  printer. 

The  Fhase  II  effort  can  be  subdivided  into  two  tasks: 

a.  C^iantify  the  brightness  of  nighttime  imagery  recorded  on  3/4”  video 
cassettes  from  an  image  intensifier. 

b.  Quantify  emitted  radiation  in  the  8-14  micrometer  spectral  band  taken 
by  a  thermal  imaging  radiometer. 
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The  Fhase  III  effort  can  be  subdivided  into  four  tasks: 


a.  Cluster  a  CTFTTAB1  data  file  and  generate  a  data  file  that  can  be  used 
as  a  background  scene. 

b.  Interactively  create  camouflage  patterns. 

c.  Overlay  patterns  onto  background  scenes. 

d.  Calibrate  the  color  monitor. 

Each  of  these  tasks  will  be  discussed  in  the  following  sections.  All  of 
the  software  has  been  implemented  and  tested  on  the  existing  Terrain  Analysis 
System. 

2.  Phase  T. 

a.  Video  Data  Processing  Software  Modifications 

The  EXPOSE  routine  creates  a  disk  file  of  reflectance  values  for  the 
wavelength  being  processed.  The  name  that  the  EXPOSE  routine  assigns  to  the 
resultant  disk  file  is  in  the  format  DLl:SCNxxx.yyy,  where  "xxx"  is  a  three 
digit  user  assigned  scene  number  and  "yyy"  is  the  wavelength  number.  The 
section  of  code  that  does  range  checking  was  modified  to  accept  for  the 
additional  values  between  720  and  900  nm. 

The  EXPOSE  routine  corrects  the  video  signal  level  of  each  pixel  within  the 
scene  area  being  processed  for  any  nonuniformities  in  the  camera  tube.  This 
corrected  value  is  then  converted  to  reflectance  by  using  the  reflectance 
target  calibration  curve  built  into  the  module.  The  reflectance  curves  of  the 
calibrated  targets  have  been  expanded  to  include  720-900  nm  data. 


1  1976  International  Commission  on  Illumination  (CIE)  L*  a*  b*  (CIELAB) . 
L*  is  a  measure  of  lightness  (0  =  black,  100  =  white) ;  a*  =  a  measure  of 
redness  (positive  axis)  or  greenness  (negative  axis) ;  b*  =  a  measure  of 
yellowness  (positive  axis)  or  blueness  (negative  axis) . 
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The  EX2DB  routine,  which  combines  the  exposure  files  into  a  data  base  file, 
has  been  expanded  to  allow  for  the  creation  of  a  data  base  file  that  contains 
27  wavelengths  in  the  range  of  380  to  900  nm. 

b.  Convert  400-900  nm  Data  to  L*  Values 

A  new  software  module  has  been  written  that  permits  the  user  to  calculate 
L*  coordinates  for  the  image  of  a  terrain  scene  as  displayed  on  near-infrared 
viewing  devices  such  as  image  converters,  image  intensifiers  and  low  light 
level  video  imaging  systems.  Ihe  user  selects  the  desired  400-900  nm  scene 
data  file  and  specifies  an  illuminant  and  the  spectral  response  of  the 
phcrtocathode  of  the  imaging  device. 

The  L*  coordinate  for  each  scene  pixel  in  the  display  is  ocnputed  in  a 
manner  analogous  to  the  calculation  of  the  visual  tristimulus  value,  Y4. 
Because  the  near- infrared  viewing  devices  present  an  image  that  varies  just  in 
lightness,  Y  and  its  corresponding  lightness  value,  L*,  are  the  only  values 
that  need  to  be  calculated.  These  values  are  calculated  in  the  following 
manner. 


A  relative  photocurrent,  i,  emitted  by  the  photocathode  of  the  imaging 
device  at  each  pixel  location  is  calculated  from  equation  (1) : 


i  = 


900 

/ 


400 


S(\)R(V)I(X)d\  anps 


(1) 


where 

V  is  the  radiation  wavelength  (micron) 

S (x)  is  the  spectral  response  of  the  photocathode  (anps/watt) 
R(x)  is  the  pixel  reflectance  spectrum  (from  the  data  base  file) 
I  ( x )  is  the  scene  illuminant  spectrum  (watts/m  -micron) 


The  Y  tristimulus  value  on  the  display  is  then  assumed  to  be  Gi,  where  G  is 
an  arbitrary  constant  that  represents  the  overall  luminous  gain  of  the  imaging 
device.  Ihe  constant  cancels  out  in  the  amputation  of  the  oolor  coordinates. 


The  tristimulus  value  Yn  of  the  white  object  stimulus  is  then  computed  from 
equation  (2) : 

Yn  "  Gin  (2) 

900 

j'  S(v)l(x)dv 

400 

This  is  the  signed,  that  would  be  generated  at  the  photocathode  by  the  image 
of  a  perfectly  diffuse  white  target.  (Reflectance  =  1.0  at  all  wavelengths.) 
The  L*  coordinate  for  each  pixel  is  computed  from  equation  (3) : 

L*  =  116/  Y  \  1/3  -16  (3) 

\M 

These  L*  values  can  then  be  clustered  to  obtain  L*  domains  for  the  displayed 
image  using  the  existing  clustering  routines. 

c.  Plotting  of  a  Single  Pixel  vs.  Wavelength 

The  software  module  REF  £C  plots  the  spectral  reflectance  curves  of  the 
pixels  closest  to  the  centroids  of  each  domain.  The  user  is  given  three 
choices  for  the  scale  of  the  reflectance  axis  -  1.0,  0.5  or  0.2.  All 
reflectance  curves  are  plotted  on  the  same  8.5"  x  11"  sheet  of  paper. 

This  software  has  been  modified  to  plot  the  additional  points  in  the  range 
720  to  900  nm. 

d.  Varying  Density  Symbol  Plots 

A  capability  has  been  provided  to  produce  varying  density  symbol  plots 
where  the  density  of  the  plotted  symbol  is  related  to  the  L*  value  of  the 
clustered  domain.  The  user  has  the  option  of  plotting  any  or  all  of  the 
domains  on  the  same  plot  or  of  separating  them  into  individual  plots  of  each 
dcmain.  These  plots,  as  well  as  all  others  produced  by  the  TAS  software,  can 


where: 

*n  = 
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be  output  on  the  Hewlett  Packard  Model  755QA  plotter  currently  attached  to  the 
TAS.  The  Hewlett  Packard  Model  7220  plotter  originally  attached  to  the  TAS 
could  also  be  used  with  no  change  in  software. 

e.  Three  Dye,  False  Color,  Infrared  Film 

In  order  to  provide  for  a  three-dimensional  clustering  process  based  on  the 
spectral  sensitivity  of  three  dye,  false  color,  infrared  (IR)  film,  a  new 
software  routine  has  been  written  that  calculates  color  values  associated  with 
each  dye  for  each  pixel  in  the  scene.  The  existing  clustering  routine,  CLUSTR, 
properly  handles  these  new  data. 

It  is  desirable  to  have  the  capability  to  evaluate  camouflage  against  the 
countermeasure  of  False  Oolor  IR  Film.  This  type  of  film  is  exemplified  by 
Kodak  Ektachrcme  IR  Film.  The  film  contains  three  layers  of  emulsion,  which 
are  sensitive  to  the  entire  visual  arri  near-IR  spectrum  to  about.  900  nm.  It  is 
normally  exposed  through  a  Wratten  Number  12  filter,  which  has  a  sharp’ lew 
wavelength  cut  off  at  500  nm. 

Upon  development,  the  exposed  silver  halide  in  each  layer  results  in  the 
appearance  of  a  proportionate  quantity  of  each  dye.  If  precisely  processed  in 
EA-5  chemicals,  the  appearance  of  the  film  will  be  proportional  to  the 
sensitivities. 

Program  FILM  reads  the  reflectance  values  from  the  data  base  image.  It 
then  calculates  the  amount  of  energy  that  would  be  incident  on  the  film  at  each 
20  nm  interval,  based  on  the  user-specified  illuminant.  Using  the  film 
sensitivity  curves,  the  program  then  calculates  the  effective  processed  dye 
value  for  each  of  the  three  dyes.  (This  is  analogous  to  the  CIE  values. ) 

These  three  values  cue  then  stored  in  a  file  for  clustering. 

f.  Domain  Plots  on  the  Digital  Letterwriter  100 

The  Digital  Letterwriter  100  is  a  dot  matrix  printer  that  is  capable  of 
producing  graphic  outputs.  A  program  has  been  written  that  displays  the 
clustered  domains  can  the  printer  as  (up  to)  5  different  densities  of  dot 
patterns.  This  produces  the  effect  of  a  grey  scale  representation  of  the 
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clustered  dona  ins.  The  user  is  able  to  specify  the  relative  dot  density 
associated  with  each  plotted  domain. 

3.  fhase_II. 

a.  Analysis  of  Image  Intensifier  Imagery 

The  TAS  is  now  capable  of  accepting  imagery  recorded  on  3/4"  video 
cassettes  taken  through  an  image  intensifier.  The  cassettes  must  be  placed  in 
the  video  cassette  recorder  (VCR)  and  played,  through  the  FOR.  A  model 

FA-410  digital  time  base  corrector.  This  is  the  identical  setup  that  is  used 
to  process  multispectral  scene  data  collected  with  the  TAS  field  data 
acquisition  system.  Additional  software  has  been  written  that  allows  the  user 
to  perform  the  desired  types  of  calculations  on  the  digitized  image  intensifier 
imagery. 

The  user  is  instructed  to  load  the  cassette  into  the  VCR  and  to  start  the 
playback  a  few  seconds  prior  to  the  beginning  of  the  frames  of  interest.  This 

gives  the  time  base  corrector  a  few  frames  in  which  to  lock  on  to  the  video 

synchronization  signed.  The  imagery  is  displayed  on  the  black  and  white 
monitor  so  that  the  user  is  able  to  see  what  frames  are  being  input  to  the 
TAS.  When  the  user  decides  that  the  proper  frame  sequence  is  being  displayed, 
a  command  to  the  software  causes  a  single  frame  to  be  digitized  and  displayed 
on  the  red/green/blue  (RGB)  monitor  (in  monochrcme) .  The  VCR  can  then  be 
stopped,  since  the  frame  of  interest  is  in  the  TAS  video  memory. 

A  digitizing  tablet  has  been  added  to  the  TAS  hardware  configuration  and 
this  tablet  is  used  to  select  image  areas  of  irregular  shape.  The  tablet  sits 

on  a  desk  next  to  the  monitor  and  has  a  "pen"  attached  to  it.  Movements  of 

this  "pen"  on  the  tablet  are  reflected  as  movements  of  a  cursor  on  the  image 
displayed  on  the  RGB  monitor.  Under  instructions  from  the  software,  the  user 
outlines  areas  of  arbitrary  shape  on  the  digitized  image  intensifier  image 
representing  a  target  and  a  background. 
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These 


A  series  of  statistics  are  provided  for  each  of  the  areas  outlined, 
statistics  include: 

Number  of  pixels  in  outlined  target 
Number  of  pixels  in  outlined  background 
Average  video  signed  level  (target) 

Average  video  signal  level  (background) 

Minimum  video  signal  level  (target) 

Minimum  video  signal  level  (background) 

Maximum  video  signed  level  (teirget) 

Maximum  video  signal  level  (background) 

Variance  of  video  signed  level  (target) 

Variance  of  video  signed  level  (background) 

Number  of  pixels  in  target  critical  dimension 
Number  of  pixels  in  height  of  video  image 
Probability  of  detection 

The  probability  of  detection,  PD,  is  a  function  of  the  signed  to  noise 
ratio  (SNR)  in  the  image5-7.  The  equations  used  in  the  software  for  these 
calculations  are  presented  in  Appendix  A. 

The  user  is  also  given  the  option  of  determining  the  relative  effect  on 
PD  that  changing  system  characteristics  produces.  The  user  is  required  to 
specify  the  conditions  under  which  data  were  collected  (CASE  1)  and  the 
proposed  new  conditions  (CASE  2) .  Data  sets  representing  these  conditions  are 
selected  from  user-maintained  data  files.  The  following  characteristics  may  be 
specified: 


Modulation  transfer  function  (MTF)  of  the  imaging  system 
Spectral  sensitivity  of  the  detector  system 
Atmospheric  transmission 
Illuminant 

Reflectance  curve  characteristic  of  target 
Reflectance  curve  characteristic  of  background 
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b.  Analysis  of  Thermal  Imagery 


Thermal  imagery  is  analyzed  in  a  manner  similar  to  image  intensifier 
imagery,  which  was  discussed  in  section  3. a.  An  additional  capability  is 
provided  in  handling  the  background.  The  user  is  given  the  option  of  either 
cxrtlining  an  area  representative  of  a  background  or  of  specifying  a  temperature 
find  emissivity.  In  addition  to  the  video  signed  level  statistics  presented  for 
image  intensifier  imagery,  statistics  are  provided  for  the  calculated 
temperature  values. 

The  provided  statistics  include: 

Number  of  pixels  in  outlined  target 
Number  of  pixels  in  outlined  background 
Average  video  signal  level  (target) 

Average  video  signal  level  (background) 

Minimum  video  signal  level  (target) 

Minimum  video  signal  level  (background) 

Maximum  video  signal  level  (target) 

Maximum  video  signed  level  (background) 

Variance  of  video  signal  level  (target) 

Variance  of  video  signed  level  (background) 

Temperature  of  average  video  signed  level  (target) 

Temperature  of  average  video  signal  level  (background) 

Tenperature  of  minimum  video  signal  level  (target) 

Temperature  of  minimum  video  signal  level  (background) 

Tenperature  of  maximum  video  signal  level  (target) 

Tenperature  of  maximum  video  signal  level  (background) 

Tenperature  equivalent  of  video  signal  level  variance  (target) 
Tenperature  equivalent  of  video  signal  level  variance  (back¬ 
ground) 

The  probability  of  detection,  PD,  is  a  function  of  the  contrast  (C) , 
contrast  threshold  (Cp)  and  the  signal  to  noise  ratio  (SNR)^-8.  The 
equations  used  by  the  software  for  these  calculations  are  presented  in  Appendix 
B. 
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The  user  is  also  given  the  option  of  determining  the  relative  effect  on 
Pq  that  changing  system  characteristics  produces.  The  user  is  required  to 
specify  the  conditions  under  which  data  were  collected  (CASE  1)  and  the 
proposed  new  conditions  (CASE  2) .  Data  sets  representing  these  conditions  are 
selected  from  user  maintained  data  files.  The  following  characteristics  may  be 
specified: 

MTF  of  the  imaging  system 

Spectral  sensitivity  of  the  detector  system 

Atmospheric  transmission 

Minimum  detectable  temperature  difference 

Emissivity  of  the  target 

Emissivity  of  the  background 

4.  Phase  III. 

The  software  for  the  Ehase  III  effort  has  been  developed  as  a  separate 
module  (DESIGN)  that  is  fully  described  in  the  User's  Manual9.  All  of  the 
capabilities  of  the  Phase  III  software  are  made  available  to  the  user  through  a 
series  of  menus. 

a.  Generation  of  Background  Scenes. 

The  Background  Scene  Generation  routine  (BAKGND)  is  used  to  create  a 
background  scene  by  clustering  a  .LAB  file.  The  routine  clusters  the 
three-dimensional  color  coordinate  data  into  a  maximum  of  250  domains.  The 
initial  clustering  is  based  upon  a  histogramming  process,  which  breaks  the 
three-dimensional  CEELAB  color  space  into  cubical  areas  that  are  tested  to  see 
if  their  pixel  density  is  greater  than  a  threshold  value.  Cubes  that  exceed 
the  threshold  cure  considered  to  be  domains.  This  process  is  repeated,  with  the 
cubes  getting  larger  on  each  pass,  until  there  are  no  more  than  250  domains. 

An  optimization  routine  may  then  be  used  to  remove  any  of  the  clustering 
inperfections  that  may  have  been  introduced  by  the  nonoptimal  histogramming 
process. 


The  background  processing  software  is  based  upon  the  clustering  software 
developed  in  previous  phases  of  this  effort.  In  these  prior  efforts,  it  was 
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desirable  to  cluster  the  scene  information  into  the  best  3,  4  or  5  domains  that 
could  be  used  for  the  camouflage  cloth.  In  this  phase,  in  order  to  generate  a 
realistic  background  scene,  it  was  desirable  to  have  as  many  domains  as 
possible,  within  the  hardware  limitations  of  the  display  system.  Since 
thehardware  display  system  inposes  a  limit  of  256  colors  displayed 
simultaneously,  and  five  colors  need  to  be  reserved  for  the  camouflage  pattern 
that  is  to  be  overlaid  on  the  background  and  one  color  (white)  reserved  for 
textual  information,  the  number  of  domains  in  the  background  scene  can  be  no 
more  than  250. 

The  BAKGND  routine  is  not  an  interactive  routine.  The  main  system  routine, 
DESIGN,  requests  information  concerning  the  .LAB  files  to  be  clustered,  as  well 
as  the  names  to  use  for  the  output  files,  and  stores  this  information  in  the 
file  BAKGND.DAT.  This  file  is  then  read  by  the  BAK3LD  routine,  which  runs  in 
the  background  and  creates  the  clustered  scene  file.  An  output  report  is 
printed  when  BAKGND  has  created  the  Background  Scene  file. 

b.  Camouflage  Pattern  Creation 

The  pattern  creation  module  (PATTON)  is  called  from  the  main  program 
(DESIGN)  when  requested  by  the  user.  The  pattern  creation  software  is  actually 
a  set  of  interactive  modules  that  may  be  selected  by  the  user  through  the  menu 
structure.  Each  of  the  individual  functions  (modules)  will  be  described  below. 

When  the  PATTON  routine  is  started  try  the  DESIGN  process,  it  initializes 
the  color  monitor  to  contain  a  blank  drawing  area  belcw  sane  descriptive 
information.  The  descriptive  information  depicts  the  color  of  each  of  the  (up 
to)  five  "pens"  that  may  be  used  to  draw  a  pattern,  an  indication  of  the  color 
of  the  current  pen,  and  a  chart  showing  the  priority  level  of  each  of  the 
colors. 

One  of  tne  first  options  that  would  normally  be  selected  is  to  specify  the 
CXELAB  values  for  each  of  the  colors  in  the  pattern  to  be  drawn.  The  user 
enters  the  number  of  the  color  to  be  set  and  its  CIELAB  coordinates.  The  RGB 
values  for  the  entered  CIELAB  coordinates  are  calculated  and  the  box  associated 
with  the  selected  color  is  filled  in  with  the  calculated  color. 
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Each  of  the  colors  that  may  be  used  to  create  the  pattern  has  a  priority 
level  associated  with  it.  This  level  is  a  user-assigned  number  between  one  and 
five,  with  each  color  having  a  unique  number.  The  level  can  be  thought  of  as 
determining  in  which  order  the  color  layers  are  applied.  First  level  one  is 
applied  to  the  "fabric",  and  then  level  two  is  applied  on  top  of  it.  This 
prooess  continues  until  all  of  the  (up  to)  five  levels  have  been  applied.  The 
priority  level,  therefore,  determines  which  color  will  be  visible  if  two  colors 
overlap  on  the  pattern.  Careful  assignment  of  level  numbers  allows  the  user  to 
draw  overlapping  patterns  without  trying  to  match  the  border  areas. 

Once  this  preliminary  information  is  entered,  the  user  may  begin  to  draw 
the  desired  pattern.  This  is  aoccnplished  by  telling  the  system  which  color  is 
to  be  applied  and  then  using  the  pen  on  the  digitizing  tablet  to  either  draw  an 
arbitrary  closed  shape  or  to  draw  a  rectangle.  As  the  pen  is  moved  on  the 
digitizing  tablet,  the  outline  of  the  shape  is  drawn  an  the  color  monitor  in 
the  selected  pen  color. 

At  any  time  during  the  pattern  creation  process,  the  user  may  select  the 
menu  option  to  fill  in  the  outlines  drawn.  The  software  uses  the  priority 
level  number  to  determine  which  color  is  visible  in  overlapping  situations. 

The  system  also  provides  an  erase  utility  that  can  be  used  to  erase 
selected  rectangular  regions  of  the  drawn  pattern. 

When  the  pattern  has  been  ccnpleted,  or  at  any  time  during  the  drawing 
process,  it  may  be  saved  to  the  disk.  In  order  to  correctly  overlay  the 
pattern  onto  a  background  scene,  however,  a  scale  for  the  pattern  Must  be 
specified.  This  is  done  by  using  the  pen  to  point  to  two  locations  on  the 
pattern  and  entering  the  distance  between  these  two  points  in  the  "fabric"  (in 
centimeters) .  The  system  also  provides  the  capability  to  save  any  portion  of 
the  drawing  area  by  allowing  the  user  to  point  to  the  upper  left  hand  comer 
and  lcwer  right  hand  comer  of  the  area  to  be  saved. 

Saved  pattern  files  may  be  recalled  for  editing.  The  full  complement  of 
pattern  creation  routines  described  above  is  available  to  edit  patterns. 
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c.  Overlaying  of  Patterns  onto  Background  Scenes. 


Patterns  created  using  this  software  may  be  overlaid  onto  clustered 
background  scenes  and  displayed  for  evaluation.  First,  the  name  of  the 
background  scene  is  entered,  and  the  scene  is  displayed  on  the  monitor.  The 
cttftar  coordinates  of  the  (up  to)  250  domains  are  translated  to  RGB  values  and 
the  background  scene  is  shewn  in  color. 

The  user  is  then  asked  to  point  to  the  upper  left  hand  and  lower  right  hard 
comers  of  the  section  of  the  background  that  is  to  be  overlaid  by  the 
pattern.  The  software  scales  the  pattern  (using  the  user-specified  scale)  to 
match  the  scale  of  the  background  scene  (determined  during  data  collection  and 
stored  as  part  of  the  file) .  The  camouflage  pattern  is  then  painted  onto  the 
screen  in  the  desired  area.  If  the  user-created  pattern  is  not  large  enough  to 
cover  the  selected  area,  the  pattern  is  repeated  in  the  horizontal  and  vertical 
directions  as  necessary. 

Onoe  the  pattern  is  overlaid  onto  the  background,  the  system  allows  the 
user  to  zoom  the  displayed  image  by  a  factor  of  two. 

d.  Monitor  Calibration 

As  described  in  the  above  sections,  background  scenes  and  overlaid  patterns 
are  displayed  in  color  on  the  monitor.  The  numeric  color  coordinates  knewn  by 
the  software  are  in  the  CIELAB  color  space.  The  monitor  needs  a  digital  number 
(0-255)  for  each  of  the  red,  green  and  blue  guns  to  display  a  color.  Software 
has  been  provided  to  transform  any  specified  L*,  a*  or  b*  coordinate  into  RGB 
values  to  produce  a  "visual  match"  to  an  object  with  the  same  L*,  a*,  b* 
coordinates.  The  fidelity  of  this  conversion  and  subsequent  display  depends 
upon  the  calibration  data  for  the  color  monitor. 

A  complete  calibration  procedure^0  is  described  in  the  TAS  Users  Manual 
and  will  not  be  repeated  here.  Appendix  C  of  this  report  presents  a  modified 
monitor  calibration  procedure.  This  short  procedure  is  not  as  rigorous  as  the 
full  calibration  procedure  and  is  useful  oily  to  display  colors  within  a  small 
area  of  CIELAB  space. 
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5.  Conclusions  and  Recommendations 


All  of  the  software  modules  for  the  three  phases  have  been  installed  and 
tested.  A  demonstration  was  held  at  Natick  using  the  Terrain  Analysis  System 
at  the  oarpletion  of  each  phase. 

A  number  of  limitations  of  the  current  system  have  become  apparent  as  more 
and  more  requirements  have  been  placed  upon  it.  It  is  therefore  reccrrmended 
that  the  following  steps  be  taken: 

a.  The  number  of  analytical  software  routines  that  form  a  part  of  the  TAS 
has  increased  dramatically  over  the  past  few  contractual  efforts.  The  slow 
speed  of  the  PDP-11/23+  processor  greatly  hampers  the  user's  efficiency  in 
performing  analyses.  The  processor  should  be  upgraded  to  allow  for  faster  and 
moire  efficient  analyses. 

b.  The  current  hardware  configuration  has  two  disk  drives,  each  of  which 
uses  removable  20  MByte  cartridges.  The  data  associated  with  each  scene  has 
increased  dramatically  as  the  system  has  been  expanded  into  the  near-IR  and 
thermal  regions.  It  is  strongly  reccnmended  that  additional  disk  storage  be 
added  to  the  system. 


This  document  reports  research  undertaken  at  the 
US  Army  Nacick  Research,  Development  and  Engineering 
Center  and  has  been  assigned  No.  NATICK, TR-  fo/otfi 
in  the  9t?ries  of  reports  approved  for  publication. 
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Appendix  A 
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Equation  set 
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Appendix  A 


PD  =  1  -  e~^SNR‘  ~ !) 
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and  I=T,B  (Target, Background) 
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The  following  terms  were  used  in  the  previous  equations: 


Pd 

SNR 

SNR' 

n*,v) 

Vr 

VB 

$ 

aB 

FACTF 

FACTM 


Nt 

Nb 

S 

T 

R 
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Nh 

Ncd 


Probability  of  Detection 
Signal  to  Noise  Ratio 
Modified  Signal  to  Noise  Ratio 

Digital  video  signal  level  for  a  pixel  at  coordinates  (x,y) 

Average  video  signal  level  for  the  outlined  target 

Average  video  signal  level  for  the  outlined  background 

Variance  of  the  outlined  target  pixels 

Variance  of  the  outlined  background  pixels 

Spectral  Sensitivity  and  Atmospheric  Factor 

Modular  Transfer  Function  (MTF)  factor 

Number  of  pixels  in  outlined  target 

Number  of  pixels  in  outlined  background 

Spectral  response  of  the  detector  (milii-amps/watt) 

Atmospheric  Transmission  (0-1) 

Reflectance  (0-1) 

IUuminant  (watts /cm7  —  micron) 

Number  of  pixel  heights  in  image 

Number  of  pixel  heights  in  critical  dimensions  of  target 


* 
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Appendix  B 

ANALYSIS  OF  THERMAL  IMAGERY 
Equation  set 
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Appendix  B 
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The  following  terms  were  used  in  the  previous  equations: 

Pd 

Probability  of  Detection 

P 

Detection  probability  of  uncluttered  targets  and  backgrounds 

SNR 

Signal  to  Noise  Ratio 

SNR’ 

Modified  Signal  to  Noise  Ratio 

r 

Ct 

Contrast  Threshold 

C 

Contrast 

c- 

Modified  contrast 

NUMCD 

Number  of  pixel  heights  in  target  critical  dimension 

f 

V1EWD 

Viewing  distance  to  display  device 

ACTWTH 

Size  of  display  device  (height) 

NUMHGT 

Number  of  pixel  heights  in  image 

V(x,y) 

Digital  video  signal  level  for  a  pixel  at  coordinates  (z,y) 

Vj 

Average  video  signal  level  for  the  outlined  target 

V'e 

Average  video  signal  level  for  the  outlined  background 

0J- 

Variance  of  the  outlined  target  pixels 

Variance  of  the  outlined  background  pixels 

FACTM 

Modular  lYansfer  Function  (MTF)  factor 

F  ACT  H 

Spectral  Sensitivity  and  Atmospheric  Factor 

FACTHC 

Spectral  Sensitivity  and  Atmospheric  Factor  (Contrast) 

FACTN 

Noise  factor 

€ 

Emissivity 

Nt 

Number  of  pixels  in  outlined  target 

Nb 

Number  of  pixels  in  outlined  background 

S 

Spectral  response  of  the  detector  (milli-amps/watt) 

T 

Atmospheric  Transmission  (0-1) 

N„ 

Number  of  pixel  heights  in  image 

Ncd 

Number  of  pixel  heights  in  critical  dimensions  of  target 

NEVD 

Noise  Equivalent  Video  Difference 

MDTD 

Minimum  Detectible  Temperature  Difference 

h 
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Appendix  C 


Calibration  of  the  DTAS  Color  Monitor 
for  CIE  L*a*b*  Color  Coordinates 

(Modified  Short  Form  Procedure) 
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Appendix  C 


This  appendix  describes  a  procedure  to  calibrate  the  DTAS 
color  monitor.  The  user  chooses  four  colors  with  known  LAB 
coordinates  that  are  representative  of  the  colors  that  are  to  be 
displayed.  Interactive  software  is  used  to  play  with  the  RGB 
values  until  a  "good"  visual  match  is  obtained  between  these  four 
displayed  colors  and  the  four  known  samples.  Off-line 
mathematical  procedures  are  then  followed  to  calculate  the 
monitor  calibration  matrix  that  can  then  be  used  by  the  software 
to  convert  an  LAB  value  into  the  appropriate  RGB  values. 

This  procedure  only  produces  "good"  matches  for  LAB 
coordinates  that  fall  within  the  area  of  the  original  four 
values.  LAB  coordinates  that  fall  outside  this  region  may 
produce  totally  false  RGB  values.  The  complete  monitor 
calibration  procedure  described  in  the  users  manual  will  produce 
visual  matches  over  a  much  greater  area  in  LAB  space. 


STEP  1. 

Select  four  color  samples  with  known  LAB  coordinates  that 
are  representative  of  the  colors  that  are  to  be  displayed.  Enter 
their  LAB  coordinates  on  lines  1-3  of  the  Linear  Calibration  Data 
Sheet  (Figure  1) . 


STEP  2. 

Using  the  monitor  calibration  software  modules  available 
through  the  DESIGN  menu  structure,  select  the  option  to  display 
uniform  boxes  of  a  user  specified  RGB  value.  Adjust  the  RGB 
values  until  the  displayed  image  is  a  "good"  visual  match  to  the 
first  known  LAB  sample.  Enter  these  values  on  lines  4-6  for 
color  one  (j=l) •  Repeat  this  process  for  the  remaining  three 
colors  and  enter  the  data  in  the  appropriate  columns  in  lines  4-6 
of  the  data  entry  form. 


STEP  3. 

Using  the  equations  shown  below  calculate  Fi ,  F2  and  F3  for 
all  four  colors  (j=l  to  4). 


Fij 


L*j+16 

116 


3 


j=l *2,3,4 
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F2j  = 


L*j+16 


116 


j=l,2,3,4 


F3j  = 


—I  3 


L*j  +16  b*j 


116 


200 


j=l,2,3,4 


Enter  these  values  on  lines  7-9  of  the  data  entry  form. 


STEP  4. 

For  each  of  the  four  colors  (j=l,2,3,4)  solve  the  following 
three  equations. 


bRx 

Fll 

bGx 

Fl2 

bBx 

= 

Fn 

Dx 

Fl4 

bRy 

F21 

bGy 

F22 

bBy 

= 

F23 

Dy 

F24 

Solve  for  bRx,  bGx,  bBx,  Dx 


Solve  for  bRy,  bGy,  bBy,  Dy 


bR2 

F31 

bGz 

F32 

bBz 

= 

F33 

Dz 

F34 

Solve  for  bRz ,  bGz,  bBz,  Dz 


i 


where: 


Vij 


Vri 

Vgi 

Vbi 

-1 

VR2 

VG2 

vB2 

-1 

VR3 

VG3 

VB3 

-1 

VR4 

VG4 

VB4 

-1 
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STEP  5. 


How  invert  the  b  matrix  as  shown  below  to  yield  U. 


Uij 


^Rx 

t>Gx 

&Bx 

*>Ry 

bGy 

bBy 

*>Rz 

bGz 

bc2 

STEP  6. 

Create  a  disk  file  using  the  editor  called  MONCAL.DAT  on 
disk  0  under  UIC  [1,15].  This  will  be  a  four  line  file  that 
contains  the  following  data: 

Line  1:  U(l,l)  U(2,l)  U(3,l) 

Line  1:  U(l,2)  U(2,2)  U(3,2) 

Line  1:  U(l,3)  U(2,3)  U(3,3) 

Line  1:  Dx  Dy  Dz 

The  format  for  the  data  is  3F10.3 
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LINEAR  CALIBRATION  DATA  SHEET 


1.  L*j 

2.  a*j 

3 .  b*j 

4.  VRj 

5.  VGj 

6.  VBj 

7.  F  !j 

8.  F2j 

9.  F3j 


Figure  1.  Data  sheet  for  linear  calibration. 


J=4 
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